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We study the exact renormalisation group flow for ultracold Fermi-gases in unitary 
regime. We introduce a pairing field to describe the formation of the Cooper pairs, 
and take a simple ansatz for the effective action. Set of approximate flow equations 
for the effective couplings including boson and fermionic fluctuations is derived. At 
some value of the running scale, the system undergoes a phase transition to a gapped 
phase. The values of the energy density, chemical potential, pairing gap and the 
corresponding proportionality constants relating the interacting and non-interacting 
Fermi gases are calculated. Standard mean field results are recovered if we omit the 
boson loops. 

The physics of ultracold fermionic gases has recently drawn much attention as it provides 
an exciting possibility of studying the regime where the dynamics of the many-body system 
becomes independent of the microscopical details of the underlying interaction between two 
fermions. This regime can be probed using the technique of Feshbach resonances [l| when 
the scattering length is tuned to be much larger then the average inter-particle separation. 
The idealised case of the infinite scattering length a is often referred to as the unitary regime 
(UR). In the a — ► — oo limit the ground state energy per particle is proportional to that of 
the non-interacting Fermi gas. 
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where M and kp are the fermion mass and Fermi momentum correspondingly and £ is the 
universal proportionality constant, which does not depend on the details of the interaction 
or type of fermions. The other dimensional characteristics of the cold Fermi-gas in the UR 
such as paring energy A or chemical potential /i can also be represented in the same way 
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H = rjE FG , A = eE FG (2) 
The infinite scattering length implies nonperturbative treatment. The most "direct" non- 
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perturbative method is based on the lattice field theory [2J, |3[ . However, being potentially 
the most powerful approach, lattice simulations still have many limitations related to finite 
size effects, discretization errors etc which may even become amplified in certain physical 
situations (the system of several fermion species is one possible example). All that makes 
the development of the analytic approaches indispensable. Several such approaches have 
been suggested so far. The (incomplete) list includes the Effective Field Theory (EFT) 
motivated formalism exploring the systematic expansion in terms of dimensionality of space 
4] and somewhat similar approach, based on 1/Af expansion Q]. More phenomenological 
approaches using the density functional method and many body variational formalism have 
been developed in {f| and 7]. 

In this paper we consider the cold Fermi gas in UR in the framework of the Exact Renor- 
malisation Group (ERG) method suggested in [Sf] and applied to the nonrelativistic many- 
fermion system with pairing in [Sj. Different aspects of the ERG approach to nonrelativistic 
systems have also been extensively studied in several subsequent papers 1(3, 11, 12, 13]. 



Being spiritually related to the EFT based approaches the ERG formalism is however fully 
nonperturbative and makes use of EFT as a guide to choose the ansatz for the effective 
action and to fix boundary conditions. The technical details of the approach in the context 



of the nonrelativistic many-fermion systems were described in 10|, ll2| so that here we give 
only a short account of the formalism. 

The central object of the ERG formalism is the average effective action which coincides 
with the bare action at the beginning of the evolution when the scale k = A (with A being 
a starting scale) and is a full quantum action when k = 0. The average effective action 
satisfies the following flow equation 



where 



d k T = - t -TT[(d k R)(T^-R)~ 1 ], (3) 



r o) = JLL-, (4) 



and Rk is the regulator satisfying the following two conditions Rk^o ~~ > and Rk-^oo ~ > k 2 . 



3 



To find a (approximate) solution of the ERG flow equations one needs to choose the ansatz 
for the effective action and fix the corresponding boundary conditions represented by the 
bare action. There are no strict quantitative criteria for choosing the ansatz so that it seems 
reasonable to work with the "relevant" degrees of freedom and include the interaction terms 
satisfying (and allowed by) all possible symmetry constraints. In our case the degrees of 
freedom are the strongly interacting fermions in the limit of the infinite scattering length. At 
the starting scale the medium parameters like Fermi-momenta, chemical potentials etc play 
a little role so that at this scale the average effective action is just the bare action with the 
standard four-fermion interaction in vacuum. One possible choice for the bare action is the 
simplest attractive EFT-motivated four-fermion pointlike interaction with the lagrangian 

A = -^C (V>W tT ) (V^)- (5) 

With a decreasing scale the role of the many-body effects becomes more and more im- 
portant and at some scale one may anticipate the Cooper instabilities, symmetry breaking, 
formation of the correlated fermion pairs etc to occur. Taking all that into account we write 
the ansatz for the effective action in the following form 

r[^t ? 0,0t] = J d * x tfw^idt + ^v^tw-u&tf) 



(6) 



We introduce the effective pairing field and used the standard Hubbard- Stratonovich 
transformation to cancel the four-fermion interaction. The couplings Z( mi M,</>,?/>) and g all 
run with the scale. We also include the kinetic term for the pairing field needed to compute 
the boson loop contributions. Note that our definition of the effective potential includes the 
term — 2fj,Z < / > (f>'(() which describes the coupling of the pairing field to the chemical potential. 
We expand the effective potential about its minimum 

U (0, f ) = u + - Po) + \ ^ 2 (0V - Po) 2 , (7) 



where the u n are defined at the minimum of the potential, <f>*<f> = po- The coefficients at the 
quadratic (in fields) term determine the phase of the system. When u± > the system is 
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in the symmetric phase with a trivial vacuum p = 0. At some critical scale the coefficient 
Ui approaches zero and the system undergoes the transition to the broken (or condensed) 
phase with p 7^ so that in this phase 

f/(0,0 t )= Mo + ^2(0 t 0-Po) 2 + , (8) 

The renormalisation factors can also be expanded about p = p 

Z<i>(<i>, f ) = z <po + %l(</>V - Po) + , (9) 

The other renormalisation factors can be expanded in the same way. 

In this paper we take into account the terms in the expansion of the effective potential 
up to quartic order in the fields. The minimum p of the effective potential evolves with 
the scale in the condensed phase and all the coefficients of the expansion should depend 
on both the scale and p . There are few options of how to organise the evolution of the 
system. The choice would depend on the physical quantities which one would like to get at 
the end of the evolution. The obvious and the most general one is just to run all quantities 
of interest. The other (simpler) way to extract essentially the same information is to fix 
some parameters like chemical potential or particle number density and run the rest. Fixed 
particle number density and evolving chemical potential seems more interesting as it gives 
the potential possibility of going to the BEC regime where the chemical potential eventually 
becomes negative. In this case the coefficients u n and the renormalisation factors will depend 
on the running scale both explicitly and implicitly via the dependencies on po{k) and p(k). 
We can, therefore, define a total derivative 

d dp d dp d 

dk =dk + lkdp-o + dkdp- (0) 

Applying this to n = —dU /dp at p = p Q and assuming the constant particle number density 
we get 

dp dp d f 

- 2 ^lk +X M=-d-p{ dk ' 

where we defined 



(11) 

P=P0 



d'u 



X= W,' (12) 

^ p=po 

The other ERG flow equations can be obtained in a similar way so that after some algebra 
we get 
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duo dp 
dk dk 



d Po . 
- M2 -dk + *** 



dU2 

dk 



U3 



dpo 
dk 



2z (f)1 



dp 
dk 

dp 
dk 



d kU\ p=P0 , 




dp 
dp 2 



d k U 



d k U 



dzA 



dk 



dp 1 , dp 
bl ~dk + 2 X dk 



1 d 2 



2 dpdp 



P=Po 



p=pa 



d k U 



P=PO 



where we have defined 



d 3 U 



X 



(13) 
(14) 

(15) 

(16) 

(17) 

dp 2 dp p=pQ 

Note that we introduce the coefficient 113 which corresponds the higher order terms in the 
expansion of the effective potential. This coefficient occurs when we act on jj^zU with the 
above defined total derivative. 

The set of evolution equations in the symmetric phase can easily be recovered using the 
fact that chemical potential does not run in symmetric phase and that po — 0. All the higher 
order terms such as XiX'i u 3 an d were calculated from the mean field type of expression 
when the boson loops are neglected and the effective potential can be calculated explicitly 
(see below). The functions po(k) and p(k) which determine the physical energy gap and 
chemical potential in the limit k — > were computed in 121 ] . In this paper we focus on the 
field independent part of the effective potential uq which is related to the energy density of 
the Fermi-gas in the UR. As can be seen from the evolution equations the coefficient uq is 
not coupled to the rest but its value depends on running chemical potential p(k) so that to 
find Uq we have to solve the whole system of the flow equations. 

The explicit expressions for the driving term d k U can be obtained from the effective 
action after the straightforward algebra and has the following form 



d k U = - 



d 3 c 



E 



FR 



(2vr) 3 



d k R F + 



27lA- 



d\ 



E 



BR 



(2tt) 3 



V 2 

B 



d k Rs- 



(18) 



where 



Ebr(q, k) = q 2 + Ui + m 2 (20V - p ) + RB(q, k), 



V B = W20V 



u 2 p 



and 



E FR (q,k,p) 



AM 

2M 



(q 2 -p 2 f) + R F(q,k,p). 



(19) 



(20) 
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The other driving term can be obtained by taking the corresponding derivatives of djJJ. The 
evolution equation as it is written above is still not enough to extract the energy density as 
it suffers from the divergence in the UV limit when k — > A. Therefore, one needs to make 
a subtraction based on a physical assumption that the energy density must be a constant, 
equal to that of the free Fermi gas when u\ = U2 = ... = u n = 0. The modified flow equation 
for the effective potential can be written as 



d 3 q E FR 1 f d 3 q E 



d k U = - / (1 ) d k R F + — / (1 " BR ) d k R B . (21) 



Note that the flow does not change as the added term is field independent. We utilise the 

n n 

type of the cut-off functions suggested first in [141] (see further discussion in [15|) for the 
boson ERG flow and in [h]] for the fermionic case. This form of the cut-off functions allows 
for the significant practical simplifications when calculating the loop diagrams. One notes 
that the cut off function for fermions can be written in different forms. However, all of them 
should contain the sgn function reflecting the particle and/or hole "faces" of the in-medium 
fermion. In this paper we use the cut off function in the form, considered in [3] 

Rf = ^ Wsgn{q - - (g 2 -pj))] 6^- \ tf-fi |), (22) 

where = (2M/i) 1//2 , and 

R B = ^(k 2 -q 2 )9(k-q). (23) 

The boson loops can be neglected at high scale and the flow equation can be integrated 
explicitly resulting in the following expression for the effective potential in the mean field 
(MF) approximation 

U MF (P, ^k) = J-0^ (E FR (p, p, k) - ^E FR{p ^ k) + fp) + C, (24) 

where C is the constant of integration. After some algebra the expression for JJ MF can be 
rewritten as 



U MF (p, p,k) = J (eM* ^ k ) + ^- q - VE FR (q,»,kr + g^ - (25) 

where e q = E FR (q, 0,0). 
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Differentiating with respect to p and setting the derivative and running scale equal to 
zero, we find that A 2 at the minimum satisfies 



M 1 f d 3 q 



4vra 2 J (2vr 



E FR (q,0,0) ^E FR (q,p F ,0y + A2 



0. (26) 



16J. 



This is exactly the gap equation derived, for example in 

To get the number density of fermions, we can differentiate U(p, /i, 0) with respect to p 
This gives 

d 3 q 



n 



(27T) 



1 _ E FR (q,p F ,0) 

v /^(g,p F ,0) 2 + A 2 



(27) 



again in agreement with Ref. [if]]. One can therefore conclude that the standard MF result 
can be reproduced within NRG if the boson fluctuations are neglected. 

The boundary conditions for the coefficients U\ can be obtained by differentiating the 



expression for U with respect to p at the starting scale k = k. 



It is worth mentioning that in this paper we include running of p, PoyU^s and the renor- 
malisation factor Z^. This is the minimum set of running parameters needed to go beyond 
the mean field approximation. The other couplings are held fixed at their initial values. 

The useful quantity to check the consistency of the approach is the boson scattering 
length clb- It is well known that the MF calculations lead to the relation = 2a F , where 
a F is the fermion scattering length. Deviation from this result is due to the boson loop effects 
and therefore goes beyond the MF approximation. In our approach the boson scattering 
length is given by the relation a# = 2uijZ\. Using the cut-off in the form specified in the 
Eqn. (22) and calculating the values of U2 and Z| in the MF approximation it can easily be 
demonstrated that the relation as = 2a F is indeed satisfied. 

The calculations with boson loops lead to the relation a# = 1.13ap. It is still quit e 
far away from the relation as = 0.6a F found in the full 4-body calculations of Ref. 17]. 
It means that the present truncation, while providing useful tool to go beyond the MF 
level, is still too crude to realistically describe the effects of boson loops at least for the 
boson scattering length. We note in passing that neither Yukawa coupling nor fermionic 
renormalisation constants run in vacuum. It is worth mentioning that similar ERG studies 
Q reSuM U1 .elation a B = 0, lo , The na tur e of the cofcence between two otherwise 
similar calculations is not clear at present. This issue clearly requires further investigations. 

The boson loop contributions are found to be small in the unitary regime when medium 
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FIG. 1: Diviation of chemical potential from its MF value as a function of pfcl. 



effects are included as the boson rescattering effects will be partially suppressed in medium. 
More definite conclusion can be drown after a number of other effects, for example evolution 
of all the couplings, is taken into account. 

Note that it is important to keep boson loops for theoretical consistency as it leads to 
the convex effective potential. The point is that in the vicinity of the physical point k = 
the coupling U2 turns out to be vanishingly small so that the potential becomes flat with 
the minimum shifted from the origin. We emphasise that the effective potential retains its 
"mexican hat" form at any finite scale. With a decreasing scale the bump becomes less and 
less pronounced so the effective potential eventually evolves into the convex form. Without 
the boson loop contribution the Ui coupling is finite at k — > and not small so that the 
convexity property of the effective potential is missing in agreement with the known results. 

In addition, smallness of the boson loop contributions in the unitary regime qualitatively 
agrees with the conclusion made in 1181 where the effects of boson loops were included via the 

n 

corresponding self-energy corrections. According to Ref. [18[ the boson loops contribution 
at zero temperature is relatively small both in the unitary regime as well as in the BCS/BEC 
regions. We show in Fig.l the deviation of the full chemical potential from its MF value as a 
function of (pfcl)" 1 in the BEC regime. This deviation is defined as £ = fi(MF) — fi(Full). 
It is seen from Fig.l this quantity is small compared to the typical value of chemical potential 
(/i ~ 0(1) at (pfcl)' 1 ~ 1). We emphasize that, although we believe that this conclusion is 
qualitatively correct, the actual value of the "beyond- mean-field" corrections can be different 
if running of all the couplings is included. Besides, including higher order terms in the 
effective action can further change the contributions from the "beyond-mean-field" effects. 
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All that constitues the subject for the future studies. We stress that, despite of being clearly 



subleading order, the boson loops contribution to the physical observables calculated in [18] 
is still larger then one obtained in this paper. To understand the cause of this deviation 
one needs to establish one-to-one diag rammatic correspondence between NRG and the more 
traditional approach adopted in 18]. It is a difficult and still unsolved problem as ERG 
includes, in principle, an infinite (and probably mixed) set of diagrams so that one-to-one 
correspondence between different contributions in the ERG approach and certain classes of 
diagrams may simply not exist. 

Now we turn to the results. First observation is that the results become practically 
independent on the starting scale A if A > 10pp. The phase transition to the condensed 
phase occurs at k cr u — Pf- One notes that the calculations with the other types of cut-offs, 
both sharp and smooth, lead to relatively close values of k cr u {9} [l2]. The fact that all values 
of k crit are clustered around the value of kp makes a good sense as at this scale the system 
becomes sensible to the medium effects like Cooper instabilities, gap formation etc. 

We found the values of 0.62,0.96 and 1.11 for the universal coefficients £,77 and e corre- 
spondingly. The calculations without boson loops give £mf(i]mf,^mf) = 0.65(0.98,1.14). 
The uncertainties, related to the form of the regulator are on the level 5%. The effect of 
the boson loops is small. It holds for both optimised cutoff function used in this paper and 
for the smoothed theta function used in Q]. One notes, however that two cutoffs lead to 
different signs of the "beyond-mean-field" contribution. In ideal case all the cutoffs should, 
of course, lead to the same results but in practise the unavoidable truncation of the effective 
action will always lead to some uncertainties. Taking into account the smallness of the bo- 
son loop effects it is hard to see if this sign uncertainty is the result of truncation or just a 
numerical instabilities which are known to be larger for expressions involving step functions. 

The obtained value for £ is close to the experimental data from (19J, £ = 0.74(7). The 
other measurements give £ = 0.34(15) and rj = 0.6(15) Q; £ = 0.32+°;^ and rj = 
2l|; £ = 0.46(5) and 77 = 0.77(5) [22J; £ = 0.51(4) and rj = 0.85(4) [2d]. The epsilon 
expansion jj] results in £ = 0.39 and 77 = 0.79. Lattice simulations give £ = 0.42(1), ^_ 



0.71(1) and e = 0.9(1) |2j, £ = 0.22(3)|24j, £ = 0.44_|25|, £ = 0.37 [26], £ = 0.3 |27|; 
£ = 0.41(2) and 77 = 0.7 [3]. Another ERG studies [l3[] give £ = 0.55. As one can see 
both experiment and numerical simulations do not provide the coherent value of the £ 
constant so it is difficult to judge the quality of the numerical estimates provided by the 
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ERG calculations. One may only conclude that the ERG approach leads to the sensible 
values of the universal coefficients consistent with the experiment and lattice calculation 
but more detailed comparison can be done when more accurate date are obtained. We note, 
however that the value of the universal parameters are still somewhat higher then the "world 



average" . One possible cause could be the neglection of the screening effects 



which are 



known to decrease the values of the gap and energy density. Naive extrapolation of our 



results using the known value of the Gorkov - Melik-Barkhudarov's correction [28| indeed 
brings the values of the universal coefficients closer to the "world average" of the lattice 
and experimental data. Clearly, this point requires further analysis and the corresponding 
work is now in progress 29] (see also jsoj]) . There are several other directions where the 
current ERG approach can further be developed. Firstlyas mentioned above, running of all 
coupling constants should be included. Secondly, the calculations with the entire effective 
potential taking into account the screening effects should also be performed 31] . Among 
the other physical applications one could mention the extension to the finite temperature 
case and analysis of the deviations from the unitary limit. 
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